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Abstract 

Charge-discharge  behavior  of  SONY  18650  lithium-ion  batteries  for  aerospace  applications  was  simulated  under  low-earth-orbit  (LEO)  condi¬ 
tions,  by  using  a  first-principles  based  mathematical  model.  The  model  determines  the  capacity  fade  on  the  basis  of  the  irreversible  loss  of  active 
lithium  ions  due  to  electrolyte  reduction.  The  capacity  fade  during  LEO  cycling  was  studied  for  5  years  of  continuous  operation  with  20%  depth 
of  discharge  as  a  function  of  the  cycling  parameters  such  as  the  end  of  charge  voltage  and  the  charging  rate. 

©  2006  Elsevier  B.V.  All  rights  reserved. 

Keywords:  Lithium-ion  battery;  Capacity  fade  simulation;  Aerospace  application;  Low-earth-orbit 


1.  Introduction 

Rechargeable  lithium-ion  batteries  have  much  higher  energy 
density,  both  volumetric  and  gravimetric,  compared  with  con¬ 
ventional  alkaline  and  nickel-hydrogen  batteries,  and  have 
been  considered  as  alternative  power  sources  to  space-qualified 
nickel-hydrogen  systems  for  orbiter  applications,  ranging  from 
‘geostationary-earth-orbit’  (GEO)  to  ‘low-earth-orbit’  (LEO) 
[!]■ 

Since  the  batteries  are  used  only  when  the  satellite  is  hidden 
from  the  sun  by  the  earth  during  eclipse  period,  the  requirements 
strongly  depend  on  the  orbit  condition  [2-5].  A  satellite  in  LEO 
(within  1000  km  from  the  earth)  uses  its  batteries  far  more  often 
than  a  GEO  satellite  (within  3.6  x  104  km  from  the  earth)  that  is 
rarely  in  the  earth’s  shadow  cone.  Thus,  the  depth  of  discharge 
(DOD)  should  be  limited  to  low  levels  of  20-40%,  in  order  not  to 
overstress  the  batteries.  Typically,  the  batteries  in  a  LEO  satellite 
require  a  life  time  over  5  years,  i.e.  the  cycle  life  exceeding  up 
to  3  x  104  cycles. 

While  lithium-ion  technology  has  been  successfully  used  in 
GEO  satellites  [6],  the  long-term  cycle  life  under  LEO  condi¬ 
tions  has  not  been  validated  yet.  There  are  only  limited  long-term 


*  Corresponding  author.  Tel.:  +1  803  111  7314;  fax:  +1  803  111  8265. 
E-mail  address:  popov@engr.sc.edu  (B.N.  Popov). 

0378-7753/$  -  see  front  matter  ©  2006  Elsevier  B.V.  All  rights  reserved, 
doi:  10. 1016/j  .jpowsour. 2006. 07. 045 


test  data  available  on  lithium-ion  batteries  under  LEO  con¬ 
ditions  [4,5].  Therefore,  model-based  estimation  of  cycle  life 
and  capacity  will  demonstrate  the  applicability  of  lithium-ion 
batteries  for  orbiter  power  sources.  The  present  work  involves 
simulation  of  charge-discharge  cycling  of  lithium-ion  batteries 
for  LEO  applications  using  a  first-principles  based  model  pre¬ 
viously  developed  by  us  at  USC  [7].  In  particular,  the  effects 
of  the  cycling  parameters  such  as  DOD,  end  of  charge  voltage 
(EOCV)  and  charging  rate  on  performance  of  the  battery  were 
investigated  under  LEO  conditions. 

2.  Simulation  model  and  cycling  protocols 

A  mathematical  model  based  on  porous  electrode  theory, 
concentrated  solution  theory.  Ohm’s  law,  and  Butler- Volmer 
intercalation  kinetics  was  developed  by  us  previously  which  pre¬ 
dicts  the  capacity  fade  of  a  lithium-ion  cell  [7,8].  In  the  model,  it 
was  assumed  that  the  capacity  fade  results  from  the  irreversible 
loss  of  active  lithium  ions  due  to  electrolyte  reduction  on 
the  anode  surface.  When  the  carbon  anode  is  polarized  more 
cathodically  than  the  reversible  potential  for  parasitic  reaction 
during  the  charging  process,  the  electrolyte  solvent  (specifically, 
ethylene  carbonate)  is  reduced  to  insoluble  salts,  resulting  in 
formation  of  a  resistive  surface  film  [9-12],  The  increase  of  the 
surface  film  resistance  with  cycling  was  also  taken  into  account 
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Nomenclature 

a\  specific  surface  area  of  porous  electrode 

(cm-1) 

CLi+  lithium  ion  concentration  in  the  electrolyte 
(mol  cm-3) 

Cp  lithium  concentration  in  the  electrode  (mol  cm-3  ) 
Cp-S  lithium  concentration  at  the  electrode/electrolyte 
interface  (mol cm-3) 

CjTiiix  maximum  lithium  concentration  in  the  electrode 
(mol  cm-3) 

D[  lithium  diffusivity  in  the  electrode  (cm2  s- 1 ) 

F  Faraday  constant  (C  mol- 1 ) 

7app  applied  current  (A) 

/,  electric  current  (A) 

Ji  total  Faradaic  current  density  across  the  electrode/ 

electrolyte  interface  (A  cm-2) 

,/p  current  density  for  intercalation/deintercalation 

reactions  (A  cm-2) 

yp  exchange  current  density  for  intercalation/ 

deintercalation  reactions  (A  cm-2) 
yp  current  density  for  parasitic  reaction  (A  cm-2) 
,/p  exchange  current  density  for  parasitic  reaction 
(A  cm-2) 

ki  rate  constant  for  intercalation/deintercalation 

reactions  (Am2,5  C-1  mol-0,5) 

<2s  loss  of  active  lithium  ions  (Ah) 

r\  radial  coordinate  (cm) 

R  gas  constant  (J  mol- 1  K- 1 ) 

Ri  radius  of  electrode  particles  (cm) 

Rs  resistance  of  newly  formed  surface  film  (Q) 

R\  total  him  resistance  (Q) 

S{  superficial  surface  area  (cm2) 

t  time  (s) 

T  absolute  temperature  (K) 

Ts  total  parasitic  reaction  time  (s) 

U  local  equilibrium  potential  (V ) 

Greek  letters 

a  transfer  coefficient  of  electrochemical  reaction 

<5  surface  him  thickness  (cm) 

s  volume  fraction  of  a  solid  phase 

4>  potential  of  a  phase  (V) 

r]  overpotential  for  electrochemical  reaction  (V) 

at  conductivity  of  newly  formed  surface  him 

(S  cm-1) 

Subscripts  or  superscripts 
1  solid  phase 

f  him  on  the  particle  surfaces 

i  negative  or  positive  electrode 

n  negative  electrode  (anode) 

N  cycle  number 

p  positive  electrode  (cathode) 

s  parasitic  reaction 


Table  1 


Key  parameters  used  in  the  simulation  model 


Parameter 

Anode 

Cathode 

Lithium  diffusivity,  D[  (cm2  s-1) 

3.8  x  to-10 

1.0  X  to-9 

Rate  constant  for  lithium 

1  x  10-s 

1  x  10-6 

intercalation/deintercalation,  k\ 

(Am2'5  C-1  mol-0'5) 

Reversible  potential  for  parasitic  reaction, 

0.38 

- 

t/socp  (V  vs.  Li+/Li) 

Exchange  current  density  for  parasitic 

2.5  x  10“12 

reaction,  /Jj  (A  cm-2) 

Conductivity  of  surface  film,  k  (S  cm-1) 

5.0  x  lO"8 

- 

in  calculation  of  the  charge-discharge  curves.  No  parasitic 
reaction  was  assumed  to  occur  on  the  cathode  surface.  The 
mathematical  model  is  fully  described  in  Appendix  A. 

In  this  work,  the  simulations  were  performed  based  on  the 
experimental  data  obtained  from  a  SONY  18650  (G8)  lithium- 
ion  cell  which  consists  of  carbon  anode  and  LiCovMnvNL02 
cathode.  The  discharge  capacity  of  the  cell  was  experimentally 
determined  to  be  2.293  Ah  at  1.126  A,  and  this  value  was  used 
for  calculation  of  charging  and  discharging  currents. 

The  previous  model  was  modified  to  simulate  the  charge- 
discharge  cycling  under  LEO  conditions  as  follows  [6]:  the 
duration  of  LEO  was  assumed  to  be  90  min,  i.e.  35  min  eclipse 
and  55  min  sunshine.  The  cell  was  galvanostatically  discharged 
for  35  min  at  a  rate  calculated  using  the  required  DOD  value. 
Next,  the  cell  was  galvanostatically  charged,  until  EOCV  was 
reached,  at  a  rate  that  is  1.1  times  the  rate  necessary  to  return 
the  discharged  capacity.  Finally,  the  cell  was  potentiostatically 
charged  at  EOCV  until  a  total  of  55  min  of  charge  was  achieved. 
The  cycle  life  was  defined  as  the  cycle  number  at  which  the  end 
of  discharge  voltage  (EODV)  dropped  below  3.0  V  [4,5].  For 
capacity  check,  the  cell  was  discharged  at  a  C/2  rate  to  3.0  V,  and 
then  it  was  charged  by  using  a  conventional  protocol  of  constant 
current  (CC)-constant  voltage  (CV),  viz.  a  C/5  rate  to  4.2  V  with 
a  50  mA  cut-off  current.  All  the  simulations  were  run  at  25  °C. 
The  key  parameters  used  in  the  model  are  listed  in  Table  1 . 

3.  Results  and  discussion 

Fig.  1(a)  shows  typical  potential  and  current  profiles  com¬ 
puted  for  the  first  LEO  cycle.  The  simulations  were  run  at 
different  DOD  values  of  20-40%  for  EOCV  of  4.05  V.  Under 
LEO  conditions,  the  orbit  duration  (charge  and  discharge  times) 
remains  constant  during  cycling,  and  thus  the  model  controls 
DOD  by  controlling  the  charging  and  discharging  currents.  It  is 
seen  that  the  larger  the  DOD  value  is,  the  lower  is  the  EODV 
value. 

Fig.  1(b)  presents  the  plots  of  EODV  of  the  cell  against  cycle 
number  measured  experimentally  and  simulated  theoretically. 
The  measurements  and  simulations  were  performed  at  20  and 
40%  DOD  values  for  EOCV  of  4.05  V.  As  a  result  of  the  irre¬ 
versible  loss  of  lithium  ions  due  to  the  parasitic  reaction  during 
the  charging  process,  the  lithium  content  (state  of  charge)  in 
the  anode  at  the  beginning  of  discharge  decreases  gradually 
with  cycling  (see  Eq.  (A.  19)).  In  order  to  maintain  the  required 
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(b)  Cycle  Number /- 


(c)  Cycle  Number  / 1 04  cycles 

Fig.  1 .  (a)  Potential  and  current  profiles  simulated  for  the  first  LEO  cycle,  (b)  the  measured  and  (c)  simulated  EODV  values  as  a  function  of  cycle  number.  The 
simulations  were  performed  with  =  2.5  x  1CT12  A  cm~2  for  EOCV  of  4.05  V. 


DOD,  the  lithium  content  in  the  anode  at  the  end  of  discharge 
decreases,  causing  the  EODV  of  the  anode  to  rise  with  cycling. 
Consequently,  the  cell  discharges  to  lower  potential,  as  the  LEO 
cycling  progresses. 

In  the  model,  the  parasitic  reaction  rate  on  the  anode  surface 
was  assumed  to  be  given  by  the  Tafel  equation  (Eq.  (A.  14)). 
The  exchange  current  density  for  parasitic  reaction  Z^1  depends 
on  the  electrolyte  composition  and  electrode  structures  in  the 
specific  battery,  and  was  used  as  an  adjustable  parameter  in  the 
simulation.  As  shown  in  Fig.  1(b),  it  was  found  that  EODV 
calculated  with  ZJ51  =  2.5  x  10-12  A  cm-2  gives  the  best  fit  to 
the  experimental  data  for  both  DOD  values. 

In  Fig.  1(c),  the  simulated  values  of  EODV  were  plotted  as  a 
function  of  cycle  number  for  EOCV  of  4.05  V  to  determine  the 
cycle  life.  From  this  figure,  the  cycle  life  was  predicted  to  be 


approximately  4.4,  2.8  and  1.9  x  104  for  the  DOD  values  of  20, 
30  and  40%,  respectively.  The  simulations  performed  using  vari¬ 
ous  cycling  parameters  indicated  that  only  the  LEO  cycling  with 
20%  DOD  in  the  EOCV  range  of  3.85^1.05  V  yields  the  cycle 
life  greater  than  3  x  104  cycles  which  meets  the  requirement  for 
LEO  applications. 

Recently,  Fellner  et  al.  [4]  reported  the  experimental  test  data 
on  the  long-term  cycling  of  lithium-ion  cells  used  for  LEO  space¬ 
crafts  up  to  9  x  103  cycles.  They  projected  the  cycle  life  of  about 
2.7  x  104  cycles  from  a  linear  extrapolation  of  EODV  toward 
3.0  V  at  40%  DOD.  However,  the  simulation  results  in  this  study 
demonstrate  that  the  plot  of  EODV  versus  cycle  number  deviates 
significantly  from  linearity,  as  the  LEO  cycling  proceeds.  The 
cycle  life  predicted  from  a  linear  extrapolation  by  Fellner  et  al. 
has  been  overestimated. 
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Fig.  2.  Galvanostatic  discharge  curves  simulated  at  different  cycle  numbers  at 
the  rate  of  C/2.  The  simulations  were  run  using  EOCV  of  4.05  V  and  DOD  of 
40%.  The  initial  discharge  curve  measured  experimentally  is  denoted  as  open 
circle  in  the  figure. 


Fig.  2  presents  the  discharge  curves  simulated  at  different 
cycle  numbers.  The  capacity  check  was  done  using  a  rate  of 
C/2.  The  LEO  cycling  was  run  at  EOCV  of  4.05  V  and  to  DOD 
of  40%.  For  comparison,  the  initial  discharge  curve  measured 
experimentally  is  also  shown  in  Fig.  2.  As  shown  in  this  figure  a 
good  agreement  between  the  experimental  and  simulated  curves 
is  observed.  The  potential  plateau  of  the  simulated  discharge 
curves  decreases  gradually  with  cycling,  which  is  attributed  to 
the  continuous  increase  in  the  resistance  of  the  surface  film.  The 
simulation  predicts  a  decrease  of  the  discharge  capacity  from 
the  initial  value  of  2.293  to  1.38  Ah  after  1  x  104  cycles. 
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Fig.  4.  Discharge  capacity  and  capacity  fade  as  a  function  of  cycle  number 
simulated  at  various  charging  rates.  The  simulations  were  performed  with  EOCV 
of  4.05  V  and  DOD  of  20%. 


Fig.  3  summarizes  the  projected  capacity  fades  of  the  cell  for 
5  years  of  continuous  operation  with  20%  DOD  and  EOCV  in  the 
range  of  3 . 85-4.05  V.  In  the  model,  the  parasitic  reaction  rate  and 
thus  the  loss  of  active  lithium  ions  depend  on  the  overpotential 
for  parasitic  reaction  )]s  on  the  anode  surface.  According  to  Eq. 
(A.  14)  in  Appendix  A,  EOCV  increases  the  overpotential  r]s 
during  the  charging  process,  thereby  accelerating  the  solvent 
reduction  reaction.  This  results  in  the  increased  capacity  fade  at 
higher  EOCV  as  presented  in  Fig.  3. 

Fig.  4  shows  the  discharge  capacities  simulated  at  different 
charging  currents.  The  simulations  were  run  using  EOCV  of 


Cycle  Number /- 


Fig.  3.  Effect  of  EOCV  on  the  capacity  fade  during  LEO  cycling  with  20%  Fig.  5.  Comparison  of  the  capacity  fades  of  lithium-ion  batteries  simulated  the- 
P)OD.  oretically  and  determined  experimentally  [6], 
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4.05  V  and  DOD  of  20%.  The  charging  currents  of  0.492,  0.541 
and  0.590  A  correspond  to  the  rates  of  0.22,  0.24  and  0.26  C, 
respectively.  As  the  charging  current  increases  from  0.492  to 
0.590  A,  the  capacity  fade  increases  due  to  the  increase  of  the 
overpotential  for  the  parasitic  reaction.  Also,  the  cycle  life  was 
predicted  to  decrease  from  45,774  to  42,373  cycles  with  increas¬ 
ing  the  charging  current,  which  suggests  that  the  use  of  low 
charging  rate  is  recommended  to  accomplish  a  long-term  for 
LEO  cycling. 

Fig.  5  compares  the  simulated  capacity  fades  with  those  mea¬ 
sured  experimentally  during  LEO  cycling  with  EOCV  of  3.95  V 
and  DOD  of  20%  [6].  The  simulated  capacity  fades  at  differ¬ 
ent  cycle  numbers  are  in  good  agreement  with  the  experimental 
values  as  shown  in  Fig.  5.  In  our  laboratory  the  experimental 
long-term  test  of  lithium-ion  cells  is  in  progress,  and  the  results 
will  be  presented  for  comparison  with  the  simulation  data  in  the 
following  paper. 

4.  Conclusions 

A  first-principles  based  mathematical  model  was  used  to 
simulate  charge-discharge  cycling  of  SONY  18650  lithium- 
ion  batteries  for  aerospace  applications  under  LEO  conditions. 
Since  the  orbit  duration  remains  constant  during  LEO  cycling, 
the  model  controls  DOD  by  controlling  the  charging  and  dis¬ 
charging  currents.  The  capacity  fade  was  calculated  on  the  basis 
of  the  irreversible  loss  of  active  lithium  ions  due  to  electrolyte 
reduction.  For  5  years  of  continuous  operation  with  a  shallow 
depth  of  discharge  of  20%,  the  capacity  fade  during  LEO  cycling 
was  studied  as  a  function  of  EOCV  and  charging  rate.  The  simu¬ 
lated  values  of  capacity  fade  were  found  to  be  in  good  agreement 
with  the  experimental  values  at  different  cycle  numbers. 
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Butler- Volmer  kinetics  was  used  to  describe  lithium  interca¬ 
lation/deintercalation 
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(A. 6) 


with 


4'  =  ¥cr  -  ct'~s) -(ct1-5) ‘(CLi+)° 
m  =  4>\  -  t/?CP  -  JiSRt 


(A. 7) 


In  Eq.  (A. 7),  it  was  assumed  that  the  variation  of  liquid  phase 
potential  along  the  current  path  is  negligible.  The  discharge 
curves  of  the  anode  and  cathode  were  measured  separately  at 
a  C/20  rate  in  a  three-electrode  electrochemical  cell  with  the 
lithium  reference  and  counter  electrodes,  and  then  the  func¬ 
tional  expressions  of  c/Pcp  were  obtained  by  the  polynomial 
regression  analysis  of  the  experimental  discharge  curves.  The 
following  condition  was  used  to  determine  the  current  decay  in 
the  potentiostatic  charging  step: 


0p  -  K  =  EOCV 

(A. 8) 

Eqs.  (A.1)-(A.4)  were  transformed  into  dimensionless  form 

9Cp  1  9  /  29Cp\ 

97  “  r?  dh  y  9?;  ) 

(A.9) 

IC  :  Cp  =  Cp,  for  0  <  r  <  R{  at  7  =  0 

(A.  10) 

9Cp 

BC  :  — —  =  0,  for  r  —  0  at  t  >  0 

9r; 

(A.  11) 
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Appendix  A 


Assuming  that  the  active  electrode  materials  are  composed 
of  uniform  spherical  particles  with  a  radius  of  //,,  the  diffusion 
of  lithium  in  the  solid  phase  is  described  by  Fick’s  second  law 
for  a  sphere 


9Cp  _  A  9  (  2dcV\ 
9 1  (r;)2  9n  y  9r;  J 


(A.l) 


The  initial  condition  (IC)  and  the  boundary  conditions  (BC) 
for  galvanostatic  charging/discharging  are  as  follows: 


IC  : 

(~<L\  _  /"fLi 

—  ^i0  ’ 

for  0  <  r; 

9Cp 

BC  : 

-  0, 

for  r[  —  0 

dn 

~  9Cp 

ip 

BC  : 

A-  '  “ 

- 1  ,  for 

9n 

F 

<  Ri  at  t  =  0  (A. 2) 

at  t  >  0  (A. 3) 

(A. 4) 


BC  : 


Cf““  A  9Cp 
— Ri 


where 


cp 


/'-'Li 

f. 

^max  ’  1 


n_ 

Ri 


Ju 

- — ,  for  r  =  1  at  t  >  0 

F 

(A.  12) 


tbi 


(A.  13) 


The  parasitic  reaction  rate  on  the  anode  surface  was  assumed 
to  be  given  by  the  Tafel  equation 

t  ,  r ,  (  oFn  F  \ 

Js  =  4  exP  y — (A-14) 

with 


bs  =  <Pl  -  A°CP  -  JnSnRl 


(A.  15) 


The  total  current  density  (7n )  at  the  anode  is  the  sum  of  the 
intercalation/deintercalation  current  density  ( /p)  and  the  para¬ 
sitic  reaction  current  density  (/p) 


1—1'- 


Li 


rLi 


for  r\  —  R,  at  t  >  0 


(A.  16) 
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The  electronic  charges  are  completely  consumed  by  the  inter- 
calation/deintercalation  of  lithium  at  the  cathode 

JP  =  (A-17) 

Op 

The  loss  of  active  lithium  ions  was  estimated  using  the  fol¬ 
lowing  equation: 

Qs  —  f  S  J^Sn  df  (A.  18) 

J  o 

Due  to  the  parasitic  reaction,  the  dimensionless  lithium  con¬ 
centration  in  the  anode  at  the  beginning  of  discharge  at  cycle 
number  (N  + 1)  is  less  than  that  at  cycle  number  N,  and  it  is 
given  by 


C; 


=  C 


Itv+i 

where 

2sL  = 


Q  s 


N 


Q 


si  JV 


SnFVnC ”« 


(A.  19) 


(A. 20) 


Continuous  precipitation  of  insoluble  product  on  the  anode 
surface  causes  the  resistance  of  the  film  to  increase  with  increas¬ 
ing  the  cycle  number 


R\ 


iV+l 


=  Ri 


R 


N 


s  I  At 


(A. 21) 


where 

R*\n=—  (A-22) 

K 
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